An SU(2) Formulation of the t-J model: Application to Underdoped Cuprates 



Patrick A. Lee 1 
Naoto Nagaosa 2 

Tai-Kai Ng 3 
Xiao-Gang Wen 1 
1 Department of Physics, MIT, Cambridge, MA 02139 
£ — ' 2 Department of Applied Physics, University of Tokyo, Tokyo 113, Japan 

3 Department of Physics, Hong Kong University of Science and Technology, Clearwater Bay, Hong Kong 
G\ ■ (June, 1995) 



We develop a slave-boson theory for the t-J model at finite doping which respect a SU(2) 
symmetry - a symmetry previously known to be important at half filling. The mean field phase 
diagram is found to be consistent with the phases observed in the cuprate superconductors, which 
contains d-wave superconductor, spin gap, strange metal, and Fermi liquid phases. The spin gap 
phase is best understood as the staggered flux phase, which is nevertheless translationally invariant 
for physical quantities. The physical electron spectral function shows small Fermi segments at low 
doping which continuously evolve into the large Fermi surface at high doping concentrations. The 
close relation between the SU{2) and the 17(1) slave-boson theory is discussed. The low energy 
effective theory for the low lying fluctuations is derived, and new lying modes (which were over 
looked in the U{1) theory) are identified. 
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I. INTRODUCTION 
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O , It is well established that high temperature superconductivity appears in cuprates when holes are doped into the 
parent compound, which is understood to be Mott-Hubbard antiferromagnetic (AF) insulators. Since the parent 
compound is insulating only by virtue of strong correlation, it stands to reason that a strongly correlated model is 
the requisite starting point to describe the cuprates. The simplest such model is the two-dimensional t-J model and 
a large effort has been made to study how the phase diagram evolves from an Heisenberg antiferromagnet when a 
concentration x of holes are introduced. The doping of a Mott-Hubbard insulator is a new problem in condensed 
matter physics and involves issues quite different from the doping of a band insulator. A key question is the evolution 
of the Fermi surface with doping. At low doping, the unit cell is doubled in the AF state and the first holes will 
\ form small pockets, not unlike the doping of band insulators. The pockets are centered on (ir/2,n/2). jjj On the 
other hand, when the hole concentration is large, it is known that a large Fermi surface is formed, with an area 
given by 1 — x, in agreement with Luttinger theorem. B The point is that the local moments on the copper are now 
counted as part of the conduction electron that makes up the Fermi sea. The key question is how this evolution takes 
place as a function of doping. It seems quite likely that the state for intermediate doping may contain novel features 
■ not encountered before. Indeed, novel concepts such as quantum spin liquid sates and spin-charge separation were 
introduced early on and much work has gone into the development of a formal theory which exhibits some of these 
novel features. j| One line of approach is to start from mean field decoupling ^J-[?J and study fluctuations about the 
mean field solution, which turns out to be an U(l) gauge theory. |§|-|Tc|l On the experimental front, much work has 
focused on the underdoped region, defined as the region of hole concentration between the onset of superconductivity 
and the maximal T c , because many anomalous properties are found in the metallic state in this regime. For example, 
unlike optimally doped systems where the magnetic susceptibility x an d the Knight shift are temperature independent, 
underdoped cuprates generally show a reduction in x below 400K or so. H] At the same time the specific heat is found 
to be suppressed relative to the T linear behavior expected for conventional metals. jl2| This behavior suggests the 
formation of a gap in the excitation spectrum. This gap also shows up in the c axis frequency dependent conductivity, 
Jk| but the conductivity in the plane is not so strongly affected. The in-plane DC conductivity shows a suppression 
below about 200K relative to the linear T resistivity observed at higher temperatures. H] This suppression can be 
attributed to reduction of the width of the Drude-like peak by a factor of two. with little effect on the spectral 
weight, the reduction of the conductivity is due to the scattering rate rather than to carrier concentration. fli"5| These 
observations suggest that the gap appears only in the spin, and not the charge degrees of freedom in the 2D plane 
and has been loosely referred to as the spin gap. We should add that the strongest gap-like behavior has been seen in 
the Cu NMR relaxation rate and in neutron scattering, both of which are sensitive to spin excitation at momentum 
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Q = (tt,tt). This latter phenomenon usually onsets at a lower temperature of order 200 K and it has been argued 
that it is observed only in bilayer or trilayer materials. jll],[L6| We shall take the point of view that the behavior at 
(71-, 7r) may be a more delicate issue depending on nesting properties at the Fermi surface etc. and for the rest of the 
paper, we shall use the term spin gap to refer to properties mentioned earlier which are characteristic of single-layer 
as well as multi-layer cuprates. 

Very recently, angle resolved photoemission experiments have yielded important new information concerning the 
electronic excitations of underdoped cuprates. It was discovered that a gap in the spectral functions already existed 
in the normal state. [|lj],|l8| Furthermore, the size of this gap and its dependence on k space is similar to the d-wave 
type gap observed in the superconducting state. The difference is that in the normal state, the gap appears to close 
in a finite segment near (tt/2, ir/2), leaving a "Fermi surface segment." If this energy gap is related to the spin gap, 
this observation gives an important boost to the notion of spin-charge separation. This is because when an electron is 
removed from the plane, as in photoemission and in c-axis conductivity, one is forced to pay the energy cost to break 
the singlet pairs in the plane, whereas for in-plane conductivity, charge transport may occur within the spin singlet 
sector. Such a behavior is in fact a natural consequence of the mean field phase diagram of the t-J model that has 
been in existence for some time. In this theory the constraint of no double occupancy is enforced by writing the 
electron operator c a i in terms of auxiliary fermions and boson particles c a i = /«,•&] and demanding that each site is 
occupied by either a fermion or a boson. In a mean field (MF) treatment, the order parameters Xij — (/<L/aj) an d 
A*j = (fufzj ~ fzifij) describes the formation of singlets envisioned in Anderson's resonating valence bond (RVB) 
picture. || Above the Bose condensation temperature of the bosons, spin charge separation occurs at the mean field 
level. In particular, in the underdoped regime the fermions are paired in a d-wave state, leading to a gap in the 
spin excitation but no gap in the charge excitation. This scenario has been used as an explanation of the spin gap 
phenomenon. [|l9 20 

While the conventional U(l) mean field theory has a lot of attractive features, it suffers from a number of defects. 
Firstly, when an attempt was made to improve the theory by including gauge fluctuations, it was found that the 
d-wave state was unstable. Q Secondly, in the underdoped regime, there are indications that the system is unstable 
to the spontaneous generation of gauge fluxes at finite wave vectors. [^2| Such instabilities will lead to a breaking 
of translation symmetry which is not observed experimentally. We note that it has recently been suggested that a 
modified d-wave state with a large gap at the (0,7r) point and vanishing gap along a segment near (tt/2, tt/2) may 
be stable against gauge fluctuations. |^| However, the question about finite wave vector instabilities remains. Such 
considerations motivated us to produce a new formulation of the constraint which generalizes the a SU (2) theory for 
half filled t-J model to t-J model away from half filling. (2^| Our hope is that since SU(2) gauge symmetry is an exact 
symmetry at half filling, the mean field approximation of the new formulation may capture more accurately the low 
energy degrees of freedom and may be a better starting point for small x. Indeed, we found that in the underdoped 
region, the mean field solution may be understood as a d-wave pairing state, or equivalently as a staggered flux 
(s-flux) phase, where the gauge flux alternates on even or odd sublattices. These states are related by local SU(2) 
gauge transformations and do not break translational symmetry. Furthermore, these states are connected smoothly 
to the 7r-flux phase at half filling which has large excitation energy at the (0, ir) point, comparable to that at (0,0). 
This is in agreement with photoemission experiments on the insulating cuprates, suggesting that the AF state may 
resemble the 7r-flux phase at short distances. |26|,^] Furthermore, in the experiment the state at (0, tt) moves towards 
the Fermi surface with doping, which can be understood in the mean field theory as a gradual closing of the spin gap. 
In this work [^4| we also introduced a residual attraction between the boson and fermions, and show that this gives 
to "Fermi surface segments" near the (tt/2, tt/2) point which grows with doping. Thus, the SU(2) mean field theory 
allows us to answer the fundamental question of how the Fermi surface evolves from hole pockets near the (tt/2, tt/2) 
point near half filling, to a large Fermi surface for large doping concentration. 

In this paper we give a more detailed description of the SU (2) theory and we also offer an alternative formulation 
which has some advantage over the original SU (2) mean field theory, particularly in the approach to large doping. 
More specifically, in the next section we show that the SU{2) theory is intimately related to the original U(l) theory. 
This leads us to a new formulation in terms of a cr-model of slowly varying boson fields. This is discussed in Sections 
3 and 4. In Section 5 we present detailed calculations of the electron spectral function, comparing the original SU (2) 
mean field approach and the new cr-model formulation. We also made some modifications of the interaction potential 
between fermions and bosons, which lead to considerable improvement of the spectral function when compared with 
experiments. In Section 6 we discuss the collective excitations of the theory, which are SU(2) gauge fields, and we 
point out the important massless gauge fields in different parts of the phase diagram. In particular, the existence 
of a massless mode in the staggered flux phase is an important new feature of the SU(2) theory compared with the 
U(l) formulation. Wc also briefly discuss the response to an electromagnetic field of the normal and superconducting 
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states. 



II. RELATION OF THE 517(2) FORMULATION TO 17(1) THEORY 

Affleck et al. p5j pointed out that the t-J model at half filling obeys an exact SU(2) symmetry. They introduced 
the SU(2) doublets 
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to represent the destruction of a spin up and spin down on site i respectively. This expresses the physical idea that a 
physical up spin can be represented by an up spin fermion or the absence of a down spin fermion once the constraint is 
imposed. The theory is invariant under the local transformation tp a i — > giip a i where g is a 2 x 2 matrix representation 
of the SU(2) group. In the original formulation, which we shall refer to as the U(l) theory, this symmetry is broken 
upon the introduction of holes. 

In Ref. pij a new formulation of the constraint of no double occupation in the t-J model was introduced which 
preserves the SU(2) symmetry even away from half filling. The key step is the introduction of a doublet of bosons, 



hi 
hi 

on each site, so that the physical electron operator can be written as an SU(2) singlet, i.e. 
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c 2i = 7/|^2; 
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The t-J Hamiltonian 
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J ( S i ■ S 3 ~ ~ A ^ n i) ~ f ( c L c «J + h - c -) 
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can now be written in terms of our fermion-boson (FB) fields. The Hilbert space of the FB system is larger than 
that of the t-J model. However, the local SU{2) singlets satisfying (IV'L'^V'cd + ^l^j)|P n y s ) — form a subspace 
that is identical to the Hilbert space of the t-J model. On a given site, there are only three states that satisfy the 
above constraint. They are f\ |0), f\\Q), and -^(b\ + ^2/2/1)!^) corresponding to a spin up and down electron, and 
a vacancy respectively. Furthermore, the FB Hamiltonian, as a SU(2) singlet operator, acts within the subspace, 
and has same matrix elements as the t-J Hamiltonian. The projection to the physical subspace is accomplished by 
introducing a set of three auxiliary fields a^, I — 1,2,3, on each site i. The partition function is written after a 
standard Hubbard-Stratonovich transformation as 



Z = J DhDht DtpDipi DaoDUe ^ 



(5) 



where the Lagrangian L is given by 



£4 [l^ai^i + hlrXj 
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The matrix 
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where Xy represents fermion hopping and Ay represents fermion pairing, respectively and J = 3J/8, £ = t/2. ]2§| ] 
The density of physical holes equals the total density of bosons 

< 1 - J m c ai >=< h\h t >=< b[bi + b\b 2 >= x (8) 

and is enforced by the chemical potential /i. 
The enforces the local constraint 

< ^L^Vm + h\r% >= (9) 

In particular, for £ = 3 we have 

< flf*i + b\ibu - blb 2i >= 1 (10) 
The Lagrangian is invariant under the local SU (2) transformation 

hi -> gjhi 
U i: j -> flftC/ijflij 

«oi^ -> ghoi^Qi - gidrgj (n) 

where ^(r) is a 2 x 2 matrix which represents an SU(2) group clement. 

Equation ^| and ^ is a faithful representation of the t-J model, just as the more standard U(l) representation is. 
The two representations must be equivalent, as long as we include all the fluctuations. To understand the relation 
between the SU(2) and the U(l) theory, we will rewrite the SU{2) theory to make it as similar to the U(l) theory 
as possible. In Appendix A we will do the reverse, i.e. we will start with the U(l) theory and write it in the form of 
the SU (2) theory, we will also discuss some subtleties of the relation. 

The key ingredient is that the two component boson field in the SU(2) representation is nothing but an SU(2) 
rotation of the standard slave boson hi, i.e. 



(12) 



The matrix gi can be parametrized as 



Zt2 %i\ 

with the constraint z* a Zi a — 1, which is satisfied by the parameterization 

zn = e 2 e 2 cos - 



9i = ( j y ) (is) 



z l2 =e l2 e l2 sin- (14) 

It is natural to introduce the iso-spin vector I 

= (sin cos <j>, sin 9 sin </>, cos 0) (15) 

Furthermore it is easy to check that 

gmgl = T-I. (16) 



4 



Thus / has the meaning of the local quantization axis parametrized by the polar coordinates 9 and <fi. The angle a in 
Zi and gi is redundant and can be absorbed into the phase of bi in Eq. [l^. Using Eq. 12 we can write Eq. [| and ^ as 



where 



Z = J DgDbDb f D^D^Da DUe fa 



V 



(17) 



l' = ~ E Tr (4^ ) + \ E ^ + «o • to«« + 

<ij> ija 



E^(^-^ + ^Tr[r 3 (. 



(?T9i) 9i) 



Ipaj 
t. 



-Tr 



(l+r 3 )gjU ijgj 



(18) 



We see that the path integral of the SU(2) theory is very similar to that of the U(l) theory. Here note that the first 
line of Eq. |l| is invariant against the local gauge transformation Eq.ll. To see this we transform of the integral 

variables ip a f = g\i) a ii 9iUijgj = Uij and gja ■ rgi — (d T gj ) gi — do • f. Then Eq. [l?] and [i"8| become 



Z = \ DgDbDtf DijiDft Da DUe 



Jo 



(19) 



and 



<i]> 



i r 3 

+ & + E 5 0^a)% + JU K 



b](d T - M + a s )bi - tJ2 Xijb]bi 



(20) 



Note that L' no longer depends on g so that the g integral can be dropped. If we drop aj' 2 integral, Eq. and |20| 
have the same form as the U(l) formulation with an exception that there t is replaced bt t = 2t. It is not our purpose 
to derive the exact equivalence between the U(l) and SU(2) path integrals, but rather we want to point out how low 
lying fluctuations in the SU(2) formulation may be reproduced in the U(l) picture. 

The U(l) mean field theory corresponds to fixing g to be unity (so that I = z) and finding and Oq ' which 

minimizes the action after summing over ip and b. In the underdoped region, it was found that U$ corresponds to 

At 



d-wave pairing of fermions. Thus the SU(2) symmetry at half filling is broken by the boson term for finite x 
the same time, it is clear that for x <C 1, there is a host of U{\) mean field states Uij = g\U^ gj which are close in 
energy to the d-wave state. Since these states are degenerate at x — 0, we may expect an energy cost of order xJ 
per hole, or x 2 J per unit cell. An example of special interest is the staggered flux phase which has a Dirac spectrum 
Ek = Vsl + A 2 , at (tt/2, 7r/2). Since the density of states of the Dirac spectrum is linear in energy, the energy cost 
is ~ /i^/AJ for a given fermion chemical potential. To satisfy the fermion number constraint, fiF ~ VxAJ so that 
in this case the energy cost is expected to be VAJ x 3 ^ 2 per unit cell. At finite temperatures, we expect that these 
low energy configurations should be included in the partition function sum. This additional degree of freedom is just 
represented by the functional integral over g in Eq. [I?], and this is the motivation for adopting the SU(2) formulation. 

In Ref. p3j a mean field theory was introduced for the SU{2) action Eq. @ and ||. The mean field is a saddle point 
of the action with respect to Uij and d , after integrating over ip, ip^ and h,fv , which is possible because the action is 
quadratic in these variables. We find that the mean field phase diagram is only slightly modified from the U(l) case, 
and consists of six different phases. 
(1) Staggered flux (s-flux) phase: 



Ui, 
Ui. 



-r 3 X + i{~) 



"A 



(21) 



and a oi = 0. In the U(l) slave-boson theory, the staggered flux phase breaks translational symmetry. Here 
the breaking of translational invariance is a gauge artifact. In fact, a site dependent SU(2) transformation 
Wi = exp[i(— \y*+ l v (7t/4)ti] maps the s-flux phase to the d-wave pairing phase of the fermions: Ui.i+x,y = — X r 3 = ' = ^' r i; 
which is explicitly translationally invariant. In the s-flux phase the fermion and boson dispersion are given by ±Ef 
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and ±Eb, where Ef — y (e/ — cig) 2 + rfj, £/ = — 2 J(cos k x +cos 77/ = — 2 J(cos fe x — cos k y )A, and a similar result 

for Eb with J replaced by t. Since iajj = we have {y ai f a i) — 1 and = (bjjjfo) = ^/2. 

(2) The 7r-flux (irF) phase is the same as the s-flux phase except here \ = A. 

(3) The uniform RVB (uRVB) phase is described by Eq. |l] with a l m = A = 0. 

(4) A localized spin (LS) phase has Uij — and a oi — 0, where the fermions cannot hop. 

(5) The d-wave superconducting (SC) phase is described by Ui i + x y = — x T 3 ± Ari and ag ^ 0, aj' 2 = 0, (61) 7^ 0, 
(b 2 ) = 0. 

(6) The Fermi liquid (FL) phase is similar to the SC phase except that there is no fermion pairing (A = 0). 

The connection with the U(l) mean field theory is now clear by using Eq. [18[ The SU(2) mean field consists of 
fixing Uij — and dg — a Q °\ For each {g,} the integral over tp,ip^, b,tf gives the free energy of a £7(1) mean field 
theory with 

u ij (g)=gtug ) g j (22) 

and 

ao • t = gja { 0) ■ fg l + g\d T g l (23) 

Upon integration over {<?i}, we see that the SU(2) mean field theory includes the U(l) mean field state {U^\a 0) } 
and all the configurations {Uij, dg} connected to it by SU(2) rotations. Thus for x <C 1 all the low energy excitations 
are included in the partition sum. This is the reason why we believe the SU(2) mean field theory is a better starting 
point for undcrdoped cuprates. 

We note that with the exception of the superconducting and Fermi liquid phases, = in the SU (2) mean field 
solution. This means that {f^Jai) — 1 and the constraint Eq. is satisfied by (b\bi) — (&I&2) = 0. Unlike the U(l) 
case, the density of fermions is not necessarily 1 — x. It is this feature which allows the staggered-flux and d-wave 
states to be gauge equivalent descriptions in the s-flux phase, for instance. One consequence is that the node in the 
gap function of the fermion excitation is pinned at (7r/2, 7t/2). In Rcf. f24f| it was found that by including an attraction 
between the boson and fermion due to the exchange of ao fluctuations, Fermi surface-like features can be recovered 
in the physical electron spectral weight which is shifted away from (tt/2,tt/2). 

A similar situation appears in the uRVB phase. The fermion Fermi surface encloses area 1 and one must go beyond 
mean field theory to produce an electron Fermi surface-like features which obey the Luttinger theorem. The problem 
is even more serious in the FL phase. Even though ag is now not equal to zero, the fermion Fermi surface area 
approaches 1 — x only very slowly with increasing x and decreasing temperature. Granted that the FL state exists 
only for x > J ft so that the motivations behind the SU(2) mean field theory is no longer applicable. Nevertheless, 
this observation means that the SU(2) mean field theory does not evolve towards the U(l) mean field theory in a way 
which is acceptable. 

We believe the origin of these difficulties lies in fixing a °^ as a mean field parameter from the beginning. For 
dg "* = 0, the constraint is satisfied on the average by (h^rh) — 0. For example, this implies (bj^i) = 0- Using Eq. [l3| 
and [l5|, this suggests that the iso-spin vector I is randomized so that (/,) = 0. On the other hand, as we approach the 
superconducting phase boundary T c from above, or the Fermi liquid boundary from the uRVB side, the boson field 
hi is becoming phase coherent and we expect that it should be slowly varying in space and time. In these regions, 
the short range correlation of the boson field is not captured by the ST/ (2) mean field theory. This motivates us to 
formulate a new effective theory for the SU (2) partition function which we shall refer to as the cr-model description. 

Our strategy is to pick a mean field configuration U^ and consider a slowly varying configuration hi in Eq. || 
or, equivalently, a slowly varying gi and b in Eq. [l8|. For each configuration, dg is solved to satisfy the constraint 
locally, after performing the integral over ip, tjy . Thus in principle dg is a functional of {hi}. Our final goal is to 
produce an effective Lagrangian for {hi} which will take the form of some nonlinear it- model to describe the low 
energy physics of the problem. This is the opposite limit to the SU(2) mean field theory: the assumption of a uniform 
dg is valid when the hi configurations are rapidly varying on the scale of the fermion correlation length, which is of 
order £g = cf/A in the s-flux phase. This picture is valid at high temperatures, whereas the cr-model approach is 
expected to be applicable near the superconducting transition and the crossover to the Fermi liquid state. The truth 
most likely lies in between the two extreme limits in most parts of the phase diagram, and it will be of interest to 
explore the consequences of both limits. 

It is clear that any tj[^ related to U^ by a 577(2) gauge transformation will give an equivalent description. Thus 
we can start with any U(l) mean field configuration. In principle, we should optimize the parameters x an( i A at the 
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end of the calculation, but in practice we expect these parameters to be not so different from that given by the U(l) 
mean field theory. We also find that a judicial choice of lljp which exhibits the symmetry of a given phase yields a 
cr-model which exhibits the proper symmetry. As a first example we discuss the uRVB state. 



III. cr-MODEL OF THE FERMI LIQUID AND THE URVB PHASES 

In U{1) mean field theory the matrix l/j^ in the uRVB state is given by llj^ = ( j5 y ^ J . Here we make the 



Xij 

choice Xij — *Xo so that ujj — ixol is proportional to the identity element. Thus itself is invariant under a 
global SU(2) transformation. 

For aj = Oq = the bosons b\ and 62 are diagonalized by the energy dispersion 

El' 2 = -2txo (sinfcz + sinfc y ) ± eft - \i (24) 

In the Fermi liquid phase, the boson condenses to the bottom of the band, located for this choice of gauge at 
Q = (7r/2,7r/2). 

As explained in Ref. J24[] the SU(2) mean field theory solution for the Fermi liquid is given by dp < 0, and 61 
contains a Bose condensed part so that (61) = bae l< ^ i ' r . Note that at finite T, thermal excitations make (^2) 7^ 0. 
From Eq. || and [l^ we see that the fermion density 

<£ /!/«*> =l-x + 2{b\b 2 ) (25) 

a 

is not equal to 1 — x so that Luttinger theorem is not obeyed. As discussed in the introduction, this motivates us to 
try the cr-model approach, where we write 

hi = hie i<3 °- p (26) 

and look for hi which is slowly varying in space and r. We can further parameterize hi = gi ^ ^ ^ . Locally we can 
consider gi = g$ as constant. By introducing ip = g\%p we see that L' in Eq. 18 takes the U(l) form 



i] a 



E 6 * ( d r - /* + J Tr fate ' f» % + |Tr [(1 + r^f ) b 3 (27) 



where 



a o ' T = 9l a o • T~gi 



(dr9 f ) 9 (28) 



The local U(l) mean field solution of Eq. |2^ is given by a' = a^z and aoo is the Fermion chemical potential chosen 
in a way which ensures that the ip fermion density is 1 — x. From Eqs. ^6] and [l^, we find that 

am = a oI(gi) (29) 
The physical electron Green's function in the SU (2) theory 

G(f,r) = -<T T (ci(r, t)c\{0,0))) = -\{T T ({h\r, t)M^ r)4(6, 0)h(6, 0))) (30) 

Assuming b is Bose condensed, we have, within the mean field theory, 

G'(r,r) = -i62(T r (/ 1 (r,r)/ 1 t (0,0))) + incoherent part (31) 

The fermion part ip a , and therefore the physical electron Fermi surface now satisfies Luttinger theorem in this slowly- 
varying approximation. 
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We would like to remark that the electron Green function in the U(l) theory has a form 

G(r,r) = -<T T (c^f,r)c\ (0,0))) = -(T T ((6^ (r, r)/ T (r, r)/|(0, 0)6(0, 0)) ) (32) 

The C/(l) mean field Green function is 

G(r,r) = -b 2 a (T T (/ T (r, r)/j (0, 0))) + incoherent part (33) 

after the boson condensation. Although the coherent part has the same dispersion relation, the quasiparticle weight 
in the U(l) mean field Green function is twice of the quasiparticle weight in the SU(2) mean field Green function. 

We next derive an explicit expression for the cr-model Lagrangian by expanding in ao and integrating out the 
fermion. This is a systematic procedure for small x. Starting from Eq. |5| the fermion integration yields a contribution 
-Trln(<9 T + JU^ } - ia ■ f). An expansion in ao to quadratic order yields the term 



£'F = ^2 ^ a o( > W n) a O /3 ( < 7^™) 7r 00 3 (°' ti; n) ( 34 ) 

where 

n$(q,uj n ) = rdre-"^e-^^^)(^ t r Q ^(r)V']^^(0)) (35) 



For the uRVB state, tt^q = Tt^Sap, where tt^ may be expanded for small q and \uj n \ < q as 

„&\q,u>) = n + C[J-V + C 2 J-*^ (36) 

The coefficient tto — —CqJ^ 1 where C , C[ and C 2 are constants of order unity. The leading term gives a contribution 
— C' J\ao\ 2 . The negative sign is a reminder that the mean field ao is a saddle point with the stable direction along 
the imaginary axis. We shall see that this negative sign yields correctly a repulsive interaction between the bosons. 
We expand in h about the bottom of the boson bands and the effective Lagrangian takes the form 

C eff = ~h)d T h + 7^-\di~ h \ 2 ~ rfJh 

+ D 1 m^ 1 {h%f + \b\ 2 a (r, r) ■ I(r, r) 

+ 2 E l a o(g,o; n )| 2 4o (o,^„) (37) 

The D\ term is used to model the repulsion between bosons and D\ is of order unity for infinite on-site repulsion. 
We have rewritten the coupling between h and a*o using Eq. |l6| and [l^. Since Eq. |37] is quadratic in ao, it can be 
eliminated, yielding a fermion contribution to the Lagrangian 

r ti ip \b\ 4 I*(q,iL> n ) ■ I(q,^ n ) 

« -^\b\ 4 J2 r *(l^n)-I(q^ n ) {-C Q J+CiJq 2 + C 2 ^\ (38) 

Using 1-1=1, the first term is iCoJ|6| 4 and it modifies the D\ term in Eq. [57|to D[ = D\ + CqJ. For J < t, this 
is a small correction. To obtain a description in terms of the z fields z = (z\,z 2 ) where Z\,z 2 are defined in Eq. [l^, 
we write h = (bo + 5b) z and integrate out the 5b field. We find 

£eff = \^Wd T z\ 2 + |6 |Vd T z + ^\d iZ \ 2 + C F (39) 
o JJ 1 Zrrib 



C f = \* 2j E r *(l^n) ■ I(q,cJ n ) (c lQ 2 + C 2 ^-\ (40) 
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We have approximated b 2 by x and &o is a constant of order -Jx at low temperature. The first term in Eq. [fli| is a 
ferromagnetic Heisenberg interaction between the isospins. Using the usual CP 1 representation, it can be written as 

^jKd.-tA^zi 2 (4i) 

where 

A i = l -{z^d i z-(d i z^z) (42) 

Note that whereas the boson part in Eq. [}9| has the full 0(4) symmetry, the fermion part only has 0(3) symmetry 
because it is independent of the overall phase a. The second term in Eq. ^ describes dissipation due to particle- hole 
excitations of the Fermi sea. Note that the fermion contribution is proportional to x 2 J which is smaller than the 
boson contribution which is proportional to xt even in the overdoped region (xt > J). For example, if T > x 2 J we can 
ignore the fermion term and if we further make the classical approximation, we conclude that at high temperatures 
the system is described by the classical 0(4) model. There is no phase transition but instead there is a cross-over 
temperature of order xt below which the phase coherence length grows exponentially. This is opposed to the U(l) 
mean field theory where there is a Kosterlitz-Thouless transition. Of course this transition is destroyed when gauge 
fluctuations are taken into account. [^9| However, in our case vortex excitations are destroyed by SU(2) fluctuations 
and we can expect a suppression of the development of phase coherence in the SU(2) formulation due to the addition 
of low lying degrees of freedom. It is interesting to ask, what is the nature of these low energy excitations? In the 
effective Lagrangian Eq. |3^, the degeneracy for constant z is a gauge symmetry: any constant z is related by a global 
gauge transformation to the U(l) uRVB state. When z is slowly varying, we can use Eq. to see that in the U(l) 
representation, Uij = ixo9i9j are generated which in general contains pairing amplitudes Ay as well as modifications 
of the hopping term \ij which affects both the boson and fermion energy. This is in contrast to the U(l) formulation, 
where only the phase fluctuation of Xij is included. Thus we may view the SU(2) formulation as a way to discover 
low lying excitations which were not so obvious in the Z7(l) picture. To complete the discussion of the low lying 
excitations we need to introduce gauge fields to the effective Lagrangian. This will be done in a later section. 



IV. cr-MODEL OF THE SUPERCONDUCTING AND THE STAGGERED FLUX PHASES 

We repeat the procedure in the last section for the staggered flux phase by choosing an appropriate J7y matrix. 

Once again any which are related by gauge transformations will give the same result, but it will be convenient 

to use a which exhibits the symmetry of the state. We have noted before that in the SU (2) mean field theory, 

the s-flux state breaks the SU(2) symmetry down to U(l). This motivates us to choose the following to describe 
the s-flux phase. 

We choose the following ansatz to describe the s-flux phase 

f#U = -*'X - HV 3 A, U i%v = + (~)^ 3 A (43) 

and 

= a o + (-)*«() (44) 
Note that U-j is invariant under global T3 rotations. In the momentum space ipi = £^ e~ l ^' ir tjjk, we have 



{ V k + a l r l , W k+Q + a l r l \ / Vfc 



where 



HLan T£, ( 4 , ^l +Q )y Wk + a^, V k+Q + a l r l J y 

H b mean = t^{hl 4+Q)(^ fe + ap i vJ+q + ) (h£ Q ) (4r,! 



Vk = -2x(sinfc 3; + sinfcy) = -2\ak 

W k = -2ir 3 A(sin k x - sin k y )= -2ir 3 Aj k (46) 
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and ^2 k ' represents summation over half of the Brillouin zone. 

To study the boson condensed phase at low temperatures, let us first assume a l = 0. In this case the boson band 
bottom is at k — (ir/2, it/2) if a l are not too large. Thus the condensed boson has a form 

h(i)\ = fh\ e -ni x+iv )„/2 ( 47 ^ 



.&a(i) J \fo 

For such a boson condensation the boson free energy is an even function of a l Q . We also note that 
satisfies 



W k +a l T L , V k+Q + a l T l J \-2iAT a -f k + a l T l , 2 X a h + a ( r' 





M(k x ,k y ,a ,a ) = \ M{k y ,k x ,a c ,-a L ) \ (49) 



Thus the fermion free energy is also an even function of a l Q . Therefore a l = is a self consistent solution. 

We would like to remark that Uij in Eq. ^3] does not contain any fermion pairing. However the boson condensate 
induces non-zero a l . A non-zero aj' 2 induces a pairing condensate of the fermions. But when aj' 2 = there is no 
pairing and the fermions are in a normal Fermi liquid state. 

Now we are ready to discuss some basic physical properties of our ansatz for different orientation of the condensate 
(&i, 62)- Without lose of the generality, we may assume &1/62 =real. In this case a 2 , = 0. We see that when b\ = 62 
(in this case <Zq = 0) the ansatz describes a translation and rotation invariant state. This state is equivalent to the 
usual ci-wave paired state in the U(l) meanficld theory after a SU(2) gauge transformation. It describes a <i-wave 
superconducting state (with a finite chemical potential) of the t-J model. When 61 7^ 62 and b\ bi 7^ 0, we have 
do 7^ and a\ 7^ 0. There is a pairing condensate in the fermions. The ansatz describes a superconducting state of 
the t-J model which also breaks the translation symmetry. The quasiparticle excitations have finite gap except at 
four isolated points near (±7r/2, ±7r/2). When &i 7^ and 62 = 0, we have 7^ and aj' 2 = 0. There is no pairing 
condensate in the fermions. The ansatz, despite the boson condensate, does not correspond to a superconducting 
state. It instead describes a Fermi liquid with broken translation symmetry and small pocket-like Fermi surfaces. 
This result is obtained through a later calculation of electromagnetic response. 

At high enough temperatures, the thermo fluctuations make (bi) = (62) = (ad) = 0. In this case the ansatz describes 
a translation and rotation invariant metallic state, which is just the s-flux phase studied in Ref. |24j 

In order to derive a tr-model for the h field, we integrate out the fermions as before. The difference now is that 
n oo = 7r oo ( W where 7Tq = 7Tq 7^ ttq . We find that 7Tq (0) « C[J whereas 7rg (0) = for a = 0. This is because 
7r z is the density-density response function and 7r z (0) is the compressibility of the fermion which vanishes due to the 
vanishing density of states in the middle of the band. For finite ajj we find that ttq O (0) — C^Oq. Now we can eliminate 
a*o to extremize the action. The problem retains rotational symmetry in the x-y plane but is anisotropic in the z 
direction. For example, for / in the x-y plane, we have Oq = and the energy of the mean field state is 

E MF = -M X x + ^x 2 (50) 

On the other hand, for I — z, we have a\ = = 0. Eliminating ajj we find the mean field energy to be 

2 7 

E MF = -4t X x + —^=x 3/2 (51) 




This result indicates that the boson condensate prefers to stay in the manifold that satisfies |6i| = I&2I, i-e., I z = 0. As 
pointed out earlier, this state is equivalent to the d-wave pairing state as opposed to the state I = z which corresponds 
to the staggered flux state with finite chemical potential. 

We can follow the procedure of the last section to derive an effective Lagrangian for the z field. The important 
difference is the appearance of the anisotropy energy. Ignoring the gradient terms from the fermion contribution, we 
can write down the following effective Lagrangian. 
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The last two terms are introduced phenomenologically to model the breaking of the 0(3) symmetry down to x-y 
symmetry when C\ ^ C3. This is adequate for / near the x-y plane but, strictly speaking, needs further modification 
near the north and south poles, due to the singular behavior of the energy cost given by Eq. 

To gain some understanding of the phases of the nonlinear cr-model, let us consider the classical limit where the r 
dependence of z are neglected. If (zi, Z2) are restricted to the manifold of minimum energy, i.e., I is in the x-y plane, 
the model is equivalent to two x-y models, with K-T transition at T^t ~ j^ nx - This temperature scale will be 

suppressed by fluctuations of / out of the x-y plane, because the energy cost per unit area is only x 2 J. However, we 
need to introduce gauge fields to Eq. [S2| before the low lying excitations can be fully discussed. 



V. THE ELECTRON SPECTRAL FUNCTION IN THE u-MODEL DESCRIPTION 



In this section we compute the physical electron Green function G(r, r) assuming that we are in the disordered 



phase of the c-model description. We have within the mean field theory 

G(r,r) = -^(T T h^(r,T)^(Q,0)h) 
w G B (r,T)G F (r,T) 



(53) 
(54) 



where 



G B (r,T) = (T T (rf(r,T)h(0,0))) 
G F (?,T) = -{T T (#{?,T)ft(p,0) 



(55) 



The boson Green function contains two parts. Note that at temperature T most bosons are in states which have 
energies of order T from the bottom of the boson band. Thus at high energies the boson Green function is given by 
the single-boson Green function G| as if no other bosons are present. The imaginary part of this part of boson Green 
function extends the whole band width of the boson band. At low energies (of order T), the boson Green function is 
determined by those nearly condensed bosons at low energies. Thus we may assume that bosons do condense and the 
second part of the boson Green function can be approximated by Const. xe l( ^ br where Qb is the momentum of the 
bottom of the boson band. From the above discussion we sec that the mean field electron Green function has a form 



Gl 0) = Const. e iQbP G F + G { °J 



The second term comes from the convolution of Gg and G F and is the incoherent part of the Green function. The 
first term is the coherent part since its imaginary part is given by discrete 5-functions. (Note those discrete (5-function 
peaks should really have a finite width of order T if the bosons do not really condense as in the s-flux and uRVB 
phases.) It is this coherent part that gives rise to the quasiparticle peaks observed in photoemission experiments. The 



(56) 



(0) 



more exact expression of G 
simplified as 



The incoherent part satisfies 



is given by Eq. |C1| in Appendix C. At low temperatures the lengthy expression can be 



Gf( w ,fc) 



(vf(k)) 2 {uf(k)Y 



EUk) 



G 



(0) 



.™ 7T A 



(57) 



(58) 



which can be shown by using Eq. CI 



In the following we go beyond the mean field theory and discuss several corrections to the mean field Green function. 
As low energies, bosons are nearly condensed. The boson fields (61, 62) (or gi(t)) change slowly in time direction (within 



11 



a range of order ^) and in spatial direction (within a range of order a-J Thus locally we may think there is really 
a boson condensation and calculate the (mean field) electron Green function in the boson condensed phase. Since in 
the different regions, the boson fields (61,62) point to different directions, the total Green function can be obtained 
by averaging the mean field Green functions for all the directions. We would like to point out that the fermion 
Green functions are different for different directions of the boson fields, because different local boson fields give rise 
to different local a l which enforces the constraint. 

The above picture of calculating electron Green function naturally comes from our er-model treatment of the SU (2) 
theory. The averaging weights for different directions are determined from the cr-model. We now make the crude 
approximation that we are in the high temperature phase of the cr-model, where all slowly varying configurations z 
are equally likely. 

We have already seen in Section 3 that this procedure yields a Fermi surface which obeys Luttinger theorem in the 
Fermi liquid phase where the bosons are condensed and in the uRVB phase where the bosons are nearly condensed. 
The fluctuations of the boson fields in the uRVB phase will give rise to finite broadening of the quasi-particle peaks. 
We now perform the same procedure in the s-flux phase. 

For each uniform configuration g t — g, ImGf (w, k) in general contains four (5-function peaks as a function of u). 
(Note for general g we have both translation symmetry breaking and fermion pairing.) After averaging over all 
orientations of g, we get a translation invariant electron Green function. This averaging also gives quasiparticle peaks 
an intrinsic width and line shape. 
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-10 1 -1 

FIG. 1. The electron spectral function for, from top down, (a) k 

k = (-tt/8, vr/8) -> (37r/8, 5tt/8), (c) k = (0, 0) -> (tt/2, tt/2), and (d) k = (0, tt/2) - 



(0,0) 



1 

-7r/4,7r/4) -> (tt/4, 3tt/4), 
We have chosen J = 1. 



(b) 



Figure |l| presents a numerical calculation of the electron spectral function using the above approximation. We have 
chosen J = J/2 and x = 1 so that the fermion band bottom is at around — 2J, to be consistent with experiments. 
We have set t = t = 2 J so that the incoherent part of the spectral function extend from — 8t = —16 J to 0, in order 
to agree with the numerical results. We have also set ^ = 0.2 so that the gap near (0,tt) is about 0.4 J. Roughly, 
the spectral function is similar to that of a <i-wave paired state with a spin gap around (0, ±7r) and (±7r, 0) of order 
A sp i„ ~ 0.4J. However the line shape and line width are quite different. If one plots ImG(w = 0, k) one can see 
that wings toward (0,tt) and (0,tt) at two sides of the peak (at (f — 6,^ — 5)) are enhanced by the averaging. It is 
because when g = X, ImG(cj = 0, k) has a Fermi pocket around (±f , ±f )■ We see that the averaging over g pushes 
the d-wave spectrum towards a spectrum which shows a segment of Fermi surface. 

In the above calculation of the spectral function, we only include some simple fluctuations (i.e. the uniform fluc- 
tuations of boson field). One may wonder how reliable is the above result? In the following we calculate ImG e by 
including some different fluctuations. We find that the spin gaps around (0, ±77) and (±tt, 0) are quite robust. However 
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the low energy spectral function near (±§, ±f ) (together with the positions and the shapes of the Fermi segments) 
are sensitive and are essentially determined by the fluctuations. Although different fluctuations have different effects, 
they in general stretch Fermi points of the mean field theory into Fermi segments. 




n y \ ✓ y 

FIG. 2. A diagram for renormalized electron Green function. The solid (dash) line is the fermion (boson) propagator. 

The dominant effect of fluctuations is to bind the bosons and the fermion into an electron. This corresponds to an 
effective attraction between the bosons and the fermions. One way to include this effect is to use the diagram in Fig. 
H to approximate the electron Green function, which corresponds to an effective short range interaction of form 



V 



with V < 0. We get 



(g1 0) )" 



(59) 



(60) 



V 



However, in general, fluctuations induce more complicated interactions. A more careful treatment can be found in 
Appendix C, where we treat two different kinds of fluctuations. The first one is the fluctuations of <Zg which induces 
the following interaction between the fermions and the bosons: 

1p^T1p ■ h)fh (61) 

The second one (whose importance was pointed out by Laughlin p0[) is the fluctuations of \Xij\ which induces 

-^ t /i) i (/i t V)i = -2tetc i (62) 

This is nothing but the original hopping term. We expect the coefficient t to be reduced due to screening, but in the 
following we adopt the form 



V(k) = U + 2£(cos k x + cos k v ) 



(63) 



for V in Eq. |6(J. Here the first and the second term comes from the first and the second kind of fluctuations. In Fig. 
H and H we plot the electron spectral function for calculated from Eq. We have chosen J = J/2, t = t = 2J, 
X = 1, A/x = 0.4, x = 0.1, and T = 0.1J. Here we choose — = 0.4 so that the renormalized gap near (0, tt) is about 
0.4J. The value U is determined from requiring the renormalized electron Green function to satisfies the sum rule 



2^ 



d 2 k 



:ImG e 



Note the mean field electron Green function in Eq. CI does not satisfy this sum rule 



J 2nJ (27T)* lmGe X/4 



(64) 



(65) 
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(71/2,71/2) (0,71) (0,0) (71/2,71/2) (71,71) (0,7C) 

FIG. 4. The points describe the dispersion of the quasi-particle peaks for the s-fiux phase in Fig. 3. The vertical bars are 
proportional to the peak values of ImGu which are proportional to the quasi-particle weight. 



We find that the gap near (0, ±7r) and (±7r, 0) survives the inclusion of gauge and \xij\ fluctuations. However 
spectral functions near (±^, arc modified. The Fermi point at (■§,§) for the mean field electron Green function 

Ge is stretched into a Fermi segment as shown in Fig. ||. We would like to point out that the electron Green function 
obtained here does not show any "shadow band" at u> = 0, i.e. ImG e (0, k) does not have any peak outside of the 
(0, 7r) - (ir, 0) line as the mirror image of the peaks that appear inside of the (0, tt) - (jr, 0) line. 
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The spectral function obtained here is qualitatively similar but quantitatively different from the one obtained in 
Ref. p4j through a similar calculation. The only difference is that here we include an additional term 2t(cosk x + 
cos k y ). Without this term the quasiparticle peaks near (0,0) get strongly renormalized and become very strong. 
The quasiparticle energies near (0, 0) get pushed so high that they are nearly degenerate with the energy gap at near 
(0,7r). Those features obviously disagree with experimental observations. After including the 2t{cosk x +cosk y ) term 
the agreement with experiments improved a lot. Due to a cancellation between the U and 2t(cosk x + cosk y ) near 
k = (0, 0), The quasiparticle energies and spectral weights near (0, 0) are quite close to the mean field values, and the 
gap at (0, 7r) now can be quite different from the quasiparticle energy at (0, 0). 

The incoherent part of the electron spectral function contains two broad peaks, each with a width about 4t. The 
incoherent part of the electron spectral function is roughly given by the boson density of states. In the SU(2) theory, 
the bosons see the staggered flux which cause the double peak structure in the boson density of states and in the 
incoherent part of the electron spectral function. As we change k, the relative weight of the two broad peaks changes 
due to the k dependence of the coherence factors u and v. The mean field results of the double-peak structure and the 
way in which the relative weight changes agree quite well with the numerical calculations. |34|] However the numerical 
calculations also observed certain shift of the positions of the two peaks as k changes. The mean field results do not 
have this shift. If we only include the U term the peak positions in the renormalized electron spectral function still 
do not shift much. However, if we include both the U and 2t(cosk x + cosk y ) terms the peak positions start to shift 
in the way similar to what observed in numerical calculations, as has been pointed by Laughlin in Ref. |30|. 



FIG. 6. A self energy diagram for the fermion Green function. The solid (dash) line is the fermion (gauge) propagator. 

The electron spectral functions calculated above have pretty sharp peaks even at high energies (say for uj w — 2J) 
in contrast to experimental findings that quasiparticle peaks are much wider at high energies. This discrepancy can 
be resolved by including the self energy of the fermion due to the gauge fluctuation. One can show that the self 
energy from the diagram in Fig. ^| is proportional to uj and k. Thus the lifetime is proportional to the inverse of the 
quasiparticle energy. To include this effect, we may assume the electron Green function to have a form 



{vf(k)Y 



( u f(k)Y 



uj - Ei(k) - ij(uj) oj-Euk) -ij(uj) 



(66) 
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If we assume the decay rate of the fermion to be j(uj) — \uj\ + 70, the resulting spectral function is quite similar to 
the line shapes observed in experiments. 

To summarize, we have considered three models which treat different types of fluctuations. 

1. Locally condensed boson picture. In this picture the quasiparticle peaks obtain intrinsic line widths and line 
shapes. Also, this picture allows us to recover the Fermi surface which has the Luttinger volume in the uRVB 
phase and the Fermi liquid phase. 

2. Short range attraction between the bosons and the fermions. Those interactions are mainly due to the Oq gauge 
and \xij\ fluctuations, which stretch the Fermi points of mean field Green function into Fermi segments. This 
attraction can also make the electron Green function to satisfy the spectral weight sum rule of the t-J model. 

3. Decay of fermions (7 oc u>) due to the gauge fluctuations. This effect broadens the quasiparticle peaks at high 
energies, and makes the spectral function look quite similar to the ones observed in experiment. 



VI. GAUGE FIELDS 



We next investigate the low lying excitations of the effective action. We first consider the uRVB state. Starting 
from Eq. B9, it is natural to introduce the transverse component of the gauge fields a,i by the standard replacement 



di — ► di + iSi ■ f + ieAi 



(68) 



Recall that in the U(l) case, transverse gauge field enforces the constraint that the sum of the fermion and boson 
current should vanish. Here the three components of the gauge field of, i — 1,2,3 enforces the vanishing of the 
analogous t £ currents corresponding to the r e density constraint given in Eq. ^[ An important difference is that in 
the SU{2) formulation the external electromagnetic field couples only to the bosons, because the physical electron 
density is given in terms of the boson density by Eq. ||, whereas in the U(l) formulation one is free to couple the A 
field to the boson or fermion, and the physical response function is the same after including the screening by the 17(1) 
gauge field, leading to the Ioffe-Larkin combination rides. We shall see how these rules are recovered or modified in 
the SU(2) case. 

In the uRVB case, it is most convenient to rotate locally to the U(l) formulation as done in Eq. ^?]. For slowly 
varying in space, we have 



U. 



= U^ + i X o(d a g}) 9j 



(69) 



where i = j + a. This is because uffl is invariant under any global rotation. The second term in Eq. |69| gives rise to 



the usual transformation property of SU(2) gauge fields: 

Oq - T = gU ■ fg + (d g^)g, 



a'-T- «t 



g'Si -Tg-i{digl)g 



(70) 



after combining with Eq. ^3[ In the rotated frame, the fermion -0 obeys the U(l) mean field solution with a chemical 
potential which enforces the fermion density to be 1 — x. We can now expand to quadratic order in a'^. The effective 
Lagrangian takes the form 



-eff 



eA ) h 



1 



ih r h + Dmi^ 1 (h^hj' 



2ml, 



-a'*(q,uj n )a" l (~q,-uj n )ir e ™(q,u) 



(fif 



ia' 4 ■ t + ieAij h 



(71) 



where h = (b, 0), 
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^™{ q ,u) = {j%,u n )j™{-q,-u n )) 



(72) 



io = ^ tr ^ 



(73) 



The spatial components are purely transverse 



(74) 



(75) 



and in the uRVB state 7r^" = Sg m Tr±. This gives rise to three degenerate transverse gauge modes which are massless 
in the uRVB state. This is confirmed by explicit calculation in appendix B. In the Fermi liquid phase b is Bosc 
condensed and these modes are massive due to the Anderson-Higgs mechanism. 

It is now clear from Eq. fn] that only the component a'^ is capable of screening the A field. Thus it is the 
component of the gauge field parallel to the quantization axis in the rotated frame which plays the role of the U(l) 
gauge field. Eq. fn] becomes 



+ 



V\b\ 2 
1 

2mi, 
1 



/(!) 



( '(2' 



9) -Un)nL(q,u) 



(76) 



The coupling of the perpendicular components £ = 1 , 2 to b may be approximated by the expectation value of 
(( a 'i ) )H 2 which simply renormalizes the chemical potential \i. The a 3 component can be integrated out and gives 
rise to the Ioffe-Larkin combination rule 



1 = (tt^) 1 + (tt£,) ] 



(77) 



where ir F is the £ = 3 component, i.e. the usual density-density response function of the fermions. Equation |7^ also 
shows that even in the SU(2) formulation, only a single component of the gauge field couples to the boson phase and 
plays an important role in suppressing the phase coherence of the boson, just as in the U(l) theory 

We next turn our attention to the s-flux phase. The main difference is that the matrix is invariant only under a 
T3 rotation gi — exp^Or^) so that the SU(2) symmetry is broken down to U(l). This produces a mass in the £ = 1,2 
modes and only the £ — 3 gauge mode remains massless. This is also checked by explicit calculation in appendix 
B. Phenomenologically we are led to the following effective Lagrangian by gauging Eq. [52| and keeping only the 
component of the gauge field. 



r 2mb \ fn |2 

C eff = —\zW z\ 



SD z - 



2nib 



\D,z\ 



4a; 2 J, 
~2Ch 



\Z\Z2 



x 2 1 



|„ |2\2 , ^ r ,Z uv "i 

\ z 2\ ) + a " 



(78) 



where 7r M1/ is the polarization tensor of the fermions for gauge field, and Dq = d T + eAa + a^r'* and Di = 
di + ieAi + ia|r 3 . 

Equation |7§| describes the low lying excitations of the underdoped regime: the superconducting and the spin gap 
phase corresponds to the ordered and thermally disordered phases of C e ff respectively. We defer a full discussion 
of this problem to a later publication. Here we give a qualitative discussion of the superfluid density in the low 
temperature phase. We integrate out a 3 in the standard way and we find the following dependence of the transverse 
electromagnetic field A*. 



C e ff = (eAi) 2 — 



{zj^zf 



1 + (mbn±/x) 



(79) 
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where 7r y = ir± (&y — -^l 2 - ) . The coefficient of A 2 is the superfluid density. On the minimum energy manifold \z± \ — \zi \ 
the second term in Eq. [79| vanishes and the superfluid density is exactly x. This is the d-wave state as we discussed 
earlier. On the other hand, when the vector / points towards the north or south pole, \z\z-i\ = and we have 
zV 3 z = ±1. In this case the fermion spectrum is that of the staggered flux phase with a finite chemical potential. 

The response n 1 *" for this orientation of / is that of a metal, which vanishes in the limit u>, q — > 0, uj < q. In this case 
the two terms in Eq. |7^ cancel and we find that p s = 0, i.e. the staggered flux phase is not a superfluid. In general, 
we expect that in the superconducting phase fluctuations of I away from the equator will cause a reduction of the 
superfluid density due to the second term of Eq. (79| For a more complete treatment, we will need more detailed 
information on -k^ v and its dependence on / which will be discussed elsewhere. 



VII. CONCLUSION 



The main result of this paper is the derivation of the effective low energy Lagrangian in terms of the boson fields 
z\ and Z2 and their coupling to gauge fields. These are given in Eq. ^ together with Eq. ^38] for the uRVB and Fermi 
liquid phases, and Eq. [78| for the s-flux and underdoped superconducting phases. In the case of the uRVB phase 
and the Fermi liquid phase, we show that the cr-model approach allows a smooth cross-over to the t/(l) mean field 
description, recovering all the desirable properties such as the Luttinger theorem for the Fermi surface area and the 
Ioffc-Larkin combination rules. This is a considerable improvement over the SU(2) mean field theory. 

In the staggered flux phase the a- model description offers some new insight into the connection of the SU{2) with 
the U(l) theory. The staggered flux phase is the disordered phase of the effective Lagrangian Eq. [7^ so that we may 
interpret the spin gap phase as fluctuations among d-wave state and s-flux states and a variety of states in between. 
While the phase diagram is quite similar to the U(l) theory, the collective excitations are very different. In the U(l) 
theory the gauge mode acquires a large mass gap of order (1 — x) J. In the SU (2) theory there are three gauge modes, 
two are massive with mass of order A, while one remains massless. We believe the low lying gauge modes may help 
stabilize this phase. In any case, the massless gauge modes will lead to large fluctuation effects which we have not 
truly explored in this paper. 

We also performed extensive numerical work to explore the consequence of the cr-model description for photoemission 
experiments. We find that within the uncertainties of the theory the qualitative features are not that different from 
the SU(2) mean field theory once the boson fermion attraction was included. We find an energy gap in the electron 
spectrum, large near (0,7r) and vanishing along a "Fermi segment" near (tt/2,tt/2). The precise size and location of 
these segments is beyond the accuracy of the present theory, but the k dependence is that of a broadened d-wave 
gap. We consider the agreement of this feature with the experiment to be strong support of the present approach. 
The SU (2) theory naturally describes an unusual superconducting transition that is not associated with opening or 
closing of spin gap. We have not treated gauge fluctuations adequately in this paper for us to describe the energy 
dependence or the lineshape of the spectral function, so that at present detailed questions which distinguishes the 
energy gap as measured from the leading edge or from the "centroid" of the spectral feature remains unanswered. 

We expect the SU (2) mean field theory to be applicable at high temperatures and the cr-model description to be 
more accurate near the phase boundary to the superconducting and the Fermi liquid phases. This is because the 
fermions respond to local fluctuations in the boson fields on a length scale of £f = J/T in the uRVB phase and 
£f = J/ A in the s-flux phase. On the other hand, the boson fluctuations are on a scale £b = {t/T) 1 ^ 2 for T > T^ E 
and £b = x- 1 / 2 for T < T { B %, where T { ° ] E = nxt is a mean field Bose condensation temperature. When £p > £g, 
we expect the fermion to average over the local boson fluctuations, and the SU (2) mean field theory is appropriate, 
whereas the cr-model approach requires that £p < £b- The difficulty is that for T < T B ° E we do not have a good 
understanding of £b because the coherence of the bosons is greatly suppressed by gauge fluctuations. In principle, we 
should solve the a model to obtain £,b to obtain a self-consistent solution, but that is beyond the scope of the present 
paper. This is why we explore the consequences of both methods and it is fortunate that the results are qualitatively 
similar in the s-flux phase. 

One important outcome of the present work is that it is clear that the transition to the superconducting state is 
very different from the conventional BCS theory. In BCS theory T c is controlled by the closing of an energy gap in the 
electronic excitation spectrum. In the present case, T c is controlled by boson fluctuations of our effective Lagrangian. 
We also note that the effective Lagrangian is not of the conventional Ginsburg-Landau form with a simple complex 
order parameter. The internal gauge degrees of freedom, parametrized by <fi and 9 (see Eq. |l4|) play an important 
role. For example, long range phase coherence can be destroyed by 9 fluctuations. Thus our picture of the normal 
phase (the disordered phase of the a model) is very different from that suggested by a number of workers, |3l] |3^] 
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based on the idea of phase fluctuations or a conventional BCS order parameter. In the latter picture normal state 
transport is due to charge 2e collective modes, whereas we have charge e metallic carriers. We believe the absence of 
signatures of strong superconducting fluctuations in the normal state favors our point of view. 

We would like to thank R.B. Laughlin for sharing his unpublished results. We also would like to thank Z.X. Shen 
for many helpful discussions. PAL is supported by NSF-MRSEC grant DMR- 94-00334 and XGW is supported by 
NSF grant DMR-94-11574 and A.P. Sloan fellowship. 



APPENDIX A: RELATION BETWEEN (7(1) AND 5(7(2) THEORY 



We start with the usual U(l) slave boson formalism where the operator c| creating an electron with spin a on site 
i is represented by the spinon (fermion) operator f} a and the holon (boson) operator 6, as 



The physical states satisfy the local constraint 

^4/ iCT + 6t6 i -l.)|phys>=Q 

Then the partition function Z of the t-J model is represented in terms of the functional integral as 

Z = J DipD^DbDb*DUDa e~~ L 



(Al) 



(A2) 



(A3) 



with the Lagrangian L being given by 



<ij> 



i r 3 



o=l 



4>ja 



L, 



(A4) 



where the first line is the Lagrangian L p for the fermions while the second line is the contribution Lb from the doped 
holes. Here the SU(2) matrix U%j has the spinon pairing order parameter Ay and the hopping order parameter Xij 

-Xli A, 



as the matrix elements, i.e., Ui- 



A* 



Xij 



The spinor il>i a is given by Eq. |l]. We have introduced three ao's. 



The 3-component ag is the time component of the U(l) gauge field corresponding to the constraint Eq. A2. The 1- 
and 2- components correspond to the constraint 



< phys'|/ii/ i2 |phys >=< phys'|//i//2lP h y s >= 



(A5) 



which are redundant [p5| and are left out in the usual U{1) formulation. 

Now we consider the SU(2) gauge transformation which is defined as the rotation of the spinor ipi in terms of a 
SU(2) matrix gi as 



ipi 

3 



Uij = gjUijgj 

a = gjaogi - {d T g\)g l 



(A6) 



The Lagrangian Lp for the fermions remains invariant with respect to the gauge transformation Eq. A6, while 
the holon contribution Lb changes. Then away from the half-filling, ([7, ao) and (U,a.o) are different configuration 
physically. Next we divide the functional integral over U and ao into two parts, i.e., the representative (U^°\ a^) 
and those which are related to it by SU(2) rotation g as g^U^g and a = gjaogi — {d T g\)gi. 
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DUDa F(U,a ) = J DU^D*^ j DgFtfU^gJ^gi-idrgbgi) 



(A7) 



No two members of (U^°\a^) are related by any 577(2) rotation g. We change the notation of the Grassmann 



variable in Eq. (A3) to ip, and then we change the Grassmann integral variables to tp = gip, ij)' — tjr to obtain 



Z = J DipDip^ DbDb* DU^ De$^ Dg exp - J drL'(ip, b, b*, C/ (0) , a { °\ g) 
The Lagrangian is given by 



(A8) 



IJOt 



(d T +a< 0) )<%+ JC/, 



(0) 

ij 



d T — n + -Tr^Ja^ - (d T g\) gi )]\ % - ^Tr[(l + T Z )g\uf 9j ] 



(A9) 



We now parameterize gi in terms of Zi using Eq. |l3| and bind the z with the slave boson b to define the SU(2) boson 
h= f£) as 



7 1Q ^2Q 



%iot bi 



This can be represented by 



hi = gi 



(A10) 



(All) 



Now the Lagrangian L' in Eq. (A9) is written in terms of b\, bi instead of b and g. First the Berry phase term is 

£ Kdrh - ^TT[r 3 (d T gl) gi }b*h 

i 

= E b i 9 T b i + ^[^2( z LdrZ ia ~ ZiadrZi^bi] 
i L a 

= ^ b*d T bi + Y z* a d T z ta b*b 
= 52b* a d T b ia (A12) 



where we have used the relation J2 a z a z a = 1 an d dr(J2a z a z a) — X) a ( z cA z ct + z^d T z a ) = 0. Next the hopping 
term of the boson is written as 



-Tr 



(1 + Ts)gluj% j 



Kb, = -thtu^hj. 



In summary the partition function Z is written as 

Z = J Dil)D^DhD}JDU {0) Da$ ) e~ So L 



with the Lagrangian L being given by 



J 



rr(0)trr(0) 
ij ij 



^E4 



•3" 



{d T +a^)5 ll+ JU\f 



i>ja 



Y i h\{d T -^ + ^)h i -tY i h\ufh j 



(A13) 



(A14) 



(A15) 
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Now the Lagrangian is invariant with respect to the SU(2) gauge transformation given in Eq. O. Then the constraint 

that no two configuration (f/W,^ ') are related by g can be relaxed because it gives only the constant gauge volume. 
Then we can drop (0) from U and a . This has the form of the SU{2) gauge theory proposed by Wen and Lee. 



However, we note that in the latter theory, the last term in Eq. A15 is replaced by t^2u ^l^ij where i = t/2. 



A possible source of this difficulty is that in Eq. A4 we impose three constraints using three Lagrangian multipliers 
ag, whereas in the standard U(l) formulation, only a single Lagrangian multiplier ao is used. We cannot justify this 
procedure because the three constraints involve noncommuting operators. Another possible source of discrepancy is 
that in going from integration over b and g to integration over h, a Jacobian may be necessary. 

APPENDIX B: MICROSCOPIC DERIVATION OF GAUGE FIELDS 

In this Appendix, we describe the microscopic derivation of the gauge fields in each of the mean field states. We 
begin by giving several arguments for when the gauge field is expected to be massless. We then show by explicit 
calculation that for the uRVB phase there are three massless transverse gauge fields. Finally we present a calculation 
of the propagator of the massless gauge field in the s-flux phase after integrating out the fermions. Because we are 
interested in the low energy dynamics, we consider only the massless gauge fields. The first task is to identify the 
massless gauge fields. For this purpose let us consider the following gauge-invariant term which appear in the free 
energy [||. 

F = Tr(Pii kUn-li r v U Vi ) (Bl) 

where 

Psj k — UijUjl ■ ■ ■ U m kUki (B2) 

is the product of U's along a closed loop i — ► j k — > i. When we write Uij as 

Uij = Ufe^ = ulfe ia v T ° (B3) 
with being the mean field configuration, we obtain the following contribution to the free energy of a"- . 



SF = TV(^L fe ^ ) e M "' T "4? t .... fc ^- w "' T °f4- ) ) 



ij k it' z f j 

= ^(U^P^^P^^P^^^ (B4) 

Then if p(°) does not commute with r a , Eq. (B4) gives the mass to the gauge field a a . For example, if P^ = C0I+C3T3 
with c's being constants, 

e' a p(°V ,a = P< 0) + i[a, P (0) ] + ^ 2 [a, [a, p(°)j] + 

= col + c 3 t 3 + C3(aV 2 - a 2 n) - y ((a 1 ) 2 + (a 2 ) 2 )r 3 + (B5) 

and a 1 and a 2 becomes massive. This is nothing but the Higgs mechanism where P is the Higgs field which are 
site variable belonging to the adjoint (vector) representation of SU(2). The condensation of P breaks the symmetry 
from SU(2) to U(l) and only one gauge field, i.e., a 3 in the above example, remains massless. On the other hand, 
if P(°) = col for every elementary plaquette, p(°) for arbitrary closed loop is const x 1 independent of the gauge 
choice. In this case we can choose a gauge where U-^ cx 1, and all the gauge fields a 1 , a 2 , a 3 remain massless. Now 
we apply the general consideration above to the each mean field state. We chose the gauge where the link variable 
for each state is given by 



uS£L = U&v=iX l (B6) 



for uRVB state, while 



(0) 

C/^ = -Xr 3 +i(-l)^A (B7) 



r(0) 
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for s-Flux and 7r-Flux states. Then the product of U's along an elementary plaquette is obtained as 



for uRVB state, while 



p(0) _ 4, 

^ P i — Xo 1 



- [(x 2 - A 2 ) 2 - 4 X 2 A 2 ]1 ± 4z X A( X 2 - A 2 )r 3 



(B8) 



(B9) 



for s-Flux and 7r-Flux states. Then it can be easily seen that in uRVB and 7r-Flux ( \ = A) states, all the gauge fields 
remain massless while only a 3 remains massless in s-Flux state. For the 7r-Flux state, we can chose the gauge where 



U 



(0) 

ii-\-x 



<(-i) <tf xi 



(BIO) 



Now we explicitly derive the effective action for the gauge fields up to the quadratic orders. We start from the 
Lagrangian in Eq. 6. We divide the link variable Uj and ao into the mean field value and the fluctuation around it. 



(0) 



a = a[, 0) + (5a 

Integrating out the fermions and bosons, we obtain the effective action for SUij and £>a . 
S cS = J Y, Tr[(t/f f + SUtHU™ + 5U i3 )] 

<ij> 

- Tr F ln{-G F0 + 5a + JSUij) + Tr B ln(-G B0 + (5a + tSUij) 



(Bll) 



(B12) 



where J = 3J/8 and Trp and Tr^ are the fermionic and bosonic traces. The Green's functions Gfo and Gbo m the 
mean field state are given by 



G F q — iu>„ 



,(°) 



G B0 = luii ~ Ms - ao - tUy 



(0) 



r(0) 



(B13) 



Now we can expand Eq. (B12) as 



S e ff = So + S\ + S*2 + 

= j E ^ 0)t ^] 



<ij> 



<ij> 



+ Tr[G F0 (<5a + JSUij)] - Tr[G S0 (^a + tSUij)] 

+ j E ^WljSUij] 



<ij> 



+ -Tr F [G F0 {5a Q + JSU tJ )G FQ (Sa + J6U V j>)] 
- ^Ti B [G B0 (5a a + tSU^G B0 {5a + t5U Vj >)] + 



(B14) 



The mean field equation is obtained from the condition that the first order terms in SU vanish. ( Here we chose the 
form Eq. (B7) for the s- and 7r-Flux states, but the obtained mean field equations are valid also for uRVB state by 
putting A = and x = Xo-) 



Jx 

JA 



d 2 k 



^2iri-2/(2jp fe )] + ^ 



4E k 



[n{-n B - 2tE k ) - n(-fi B + 2tE k )] 



[1 - 2f(2JE k )} + 



tAri 2 

—J±[n(-n B - 2tE k ) - n{-n B + 2tE k )} 



(B15) 
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where 



7fc = cos k x + cos k y 
rjk = cos k x — cos k v 



(B16) 



and /(&) and n(x) are the Fermi and Bose distribution functions, respectively. The condition that the first order 
terms in Sao vanish gives Eq. ^, which is satisfied by the mean field solutions. 

Now we study the second order terms S^- From the considerations given above we consider only the gauge fields 
which commute with . First consider the uRVB and 7r-Flux states, where oc 1 and all the gauge fields remain 
massless. When we make a gauge transformation 



(o) 



U, 



the action does not change. Let us take g t = e l 



If gi commutes with U-j , we have 



[7(0)^(0,-^)7 



(B17) 



(B18) 



If we consider this as an expansion SUij about Jj[^ which corresponds to a pure gauge configuration, we can see that 
the coefficient of the second order term in a° = — Oj vanish for any a = 1,2,3 for the uRVB and 7r-Flux phases, 
and only for a = 3 for the s-Flux phase. 

Generally the second order contribution 52 can be written as 



S * = J EE E n*(q,iu n )al(q,iu; n )aZ(-q,-iu; n ) 



a.b fl,u q,iuj n 

and the above consideration guarantees the masslessness of a a , and leads to the condition 

n°t(g = o,w„ = o) = o. 



(B19) 



(B20) 



Here a,b = 1,2,3 for the uRVB and tt-FIux states while a = b = 3 for the s-Flux state. 

To make things more clear, we describe here the explicit calculation for the uRVB state. Taking the gauge choice 
of Eq. (B6), the mean field equation is obtained as 



1 



d 2 k 



2 (2tt)2 



(B21) 



where 7^. = sink x + sink y . At first glance this appears different from Eq. (B15), but it can be shown by using partial 
integration that Eq. (B15) is reduced to Eq. (B21) by putting A = and \ = \a- This can be also written as 



Jxo = - 



d 2 k 



J'Xo%f'{-2Jxolk) + t 2 xvT k n'{-VB - 2t X olk) 



(B22) 



where f'(x) = df(x)/dx (n'(x) — dn(x)/dx ). 

Now we consider the second order contribution $2- The gauge fields are related to SU^ as 

8U ll+ll = 6U} i+fl = -xoof i+At T 
Then the coupling 5; n t between the fermions (bosons) with the gauge field is written as 



S' u 



a, \iv kq 



COS k,. 



2Jxoa'] 1 (q)il>k+q/2Taipk-q/2 + 2txoa a ^(q) hk+q/2 T ahk-q/7 



where ipi^ are spinors. Now S2 is explicitly given by 

S * = E Yfi J A» ~ K» b (l) - n*?(q)}a°(q)a b v (-q) 



(B23) 



(B24) 



(B25) 



a, [iv q 
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where 

Kf(l)=^»J 2 X 2 o f ^cosk.cosk j^- 1 ^-^ (B26) 

J--K ( 27r ) tuJ rn - Kk+q/2 + Kk-q/2 

where £fc = — 2JxoJk- Similar expression is obtained for Il^ b (q). It can be easily seen that in the limit q — ► 0, 
n>" fc (<7) + n^" fc ((7) — > 2 Jx^S^v by using the mean field equation Eq. (B22). The similar cancellation is obtained for 
a 1 , a 2 , a 3 in the 7r-Flux state, and for a 3 in the s-Flux state. 

Finally we present a calculation of the gauge field propagator when the fermions are integrated out, it we compute 
Il^ b (q,ioj n ) — H^ b (q,iuj n ) — U^ b (q = 0,ioj n = 0) in terms of the continuum approximation in the limit of small 
vpq,uj n , and T compared with J. For the uRVB state, this calculation is exactly the same as in the U(l) case 
described in ref.10. For the s-Flux state, we consider the following effective Lagrangian for the fermions in the 
continuum approximation: 

L = J d 2 ryj\[D T - 2iJ X (D x + D y )a 3 r 3 + i2JA(D x - A> 2 ]V>i 

+ J d 2 r^\ [D T - 2iJ X (D x - D y )a 3 T 3 + i2JA(D x + A> 2 ]V>2 (B27) 

where = + ifi 3 . Since the problem has relativistic symmetry, it is convenient to introduce a M = (— ia®, a x , a y ). 
The spinor ipi (^2) describes the fermions near ±(7r/2, 7r/2) (±(7r/2, — 7r/2)) in the fc-space. Because of the double 
periodicity in the gauge choice of Eq. (B7), k and k+ (ir, n) are coupled and cr's are the Pauli matrices describing this 
2x2 space in addition to the original isospin space spanned by r matrices. When we integrated over the fermions in 
Eq. (B27), the following integral g M1 /(q) appears. 

^ (q) -J(^kfTf (B28) 

where k = (fc , ki,k 2 ) = (u, k) is the vector in (2+l)-dimensions. 
By using the Feynman's trick, i.e, 

1 f 1 dz .„ 

(B29) 



the integral is transformed as 



ab J a [a? + 6(1- z)} 2 ' 

d 3 k fc M (k + q)„ 



(2tt)3 [(k + zq) 2 + z(l-z)q 2 ] 2 



I d*k k 6 + z(l-z)q q„ 

J (2tt) 3 [k 2 + z(l - z)g 2 ] 2 V ' 

Now (q) is diverging if the ultraviolet cut-off A for fc-integration is infinity. This is cured if one consider (q) — 
.9^(0), which is converging when A — > 00. Using J A^ WT^+^ = 5^ where q = q 2 + u 2 , we obtain 

9,M) ~ = ^ =^ [ dzjzjl^z) = (B31) 
Using Eqs. (B31) and (B32) we obtain the following effective action of the gauge field at zero temperature. 



5 e ff = y~ 



2 _l A2 

(B32) 



J 2 X A/ xy (q)/ xl/ (-q) + X 16 ^f (/o*(q)/o*(-q) + /o^/o^-q)) 



where / M „ = 9 M a„ — <9„a M . The coefficient of a^a v of this expression gives the inverse of the gauge propagator which 
is correct for small q and uo in the lattice model. 
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APPENDIX C: ELECTRON SPECTRAL FUNCTION 



The more exact expression of mean field electron Green function G^ is given by 
G<°> ee (T (c T c{)) = 11 ]T - + u ^ - *)"'(«)] 2 

x E K(^(5-fc))+n F (^(g))] 



<*=+,- 
1 1 



u — 



El(q)-E* a {q-k) 



-iS 



<*=+, 



Here N is the number of sites and 



x Y [n b {E b a (q-k))+n F (E f _ a (q)) 



El{q)-E h _ a {q-k) 



iS 



(CI) 



Ei{k)=±^(e{) +(4) 

e{ = 2 Jx(cos k x + cos k y ) 
tjI = 2JA(cosfc r — cosfc y ) 



E b ± (k)=±^(eiy + (4Y-f, B 
e\ = 2tx(cosk x + cos k y ) 
rjl = 2tA(cosk x — cos k v ) 



u f < b (k) = 



v f ' b (k) = 



1 

V2' 



J* 



1 vi 



f,b 



ribj(E) are boson or fermion occupation numbers at energy E. The incoherent part comes from the terms containing 
rif (E±). One can show that ImG-° = for lu > and 



— ImG<°> = - 
2ir m 2 



(C2) 



The coherent parts come from the terms containing nb{E h _). Note rib (E b _{k)) is almost zero except near (0,0) and 
(it, 7r). Approximating those peaks by S- functions in /c-space, we get 



(yf(k)) 2 | (^(fc)) 2 
oj-E f _(k) Lu-El(k) 



(C3) 



Next we will consider effects of fluctuations. We will consider two different types of fluctuations. The first one is 
the fluctuations of whose effect is modeled by the following effective short range interaction between the fermions 
and the bosons: 



tp^fiph^rh 

The second one is the fluctuations of \xij\ which induces 



-2tac] 



(C4) 



(C5) 



In the s-flux phase the electron operator c 1 — -j=tp^h mixes with an operator c) = -^ip^r 3 h. We find (Tc^qcI) 



iG*q with Q = (tt, it) is nonzero and given by 
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G Q {k) = ^ ~ k ) + v f (q)v b (q - k)] x [vf(q)u b (q - k) - u^q)v\q - k)\ 

1 

£ <*/3[n b (E%(q-k))+n F (Efa))]- 

a,/3=+,- <*■ 



EL{q)~E b Jq-k) 



iS 



(C6) 



—v^\k)u^(k) 



u-EUk) uj-E f _(k) 



G% 



(C7) 



Since the interaction couples to both c and c, we have to invert a two-by-two matrix to calculate the electron Green 
function. Noticing that (cfeCfe) = {c\cy and introducing 



G = 



I Gf\k) -iG$\k) 



G^ik) G { °\k + Q) 



we find that the electron Green function is the (1,1) component of the two- by- two matrix 

n -1 



G e 



G 



U x 

u 2 



(C8) 



(C9) 



Note that when Ui = the above equation reduces to Eq. |6C|. XJ\p are obtained by rewriting the interaction Eq. C4 
in the c and c basis: 

U^fiph^fh = WJc - U^rhtffip (CIO) 
= WJc-U$c+ ... (Cll) 
term has a form ip^T^-^hh^r 1 ' 2 -^ and does not contribute to the electron Green function. Thus U\ = 3/7 



-U for the interaction in Eq. C4 



The 

and U2 = 

At low energies the interactions are dressed by fermion bubbles. Because ajj gauge field is massless in the s-flux 
phase, its fluctuations mediate a long range interaction. Thus the interaction ip^T 3 ij}h''T 3 h is enhanced at low energies. 
To study this effect let us consider an extreme case which has the following interaction 



Uc^c + C/c t c + ... 



(C12) 
(C13) 



The term does not contribute to the electron Green function. Thus the electron Green function in this case is 
given by Eq. [c| with Ui = U 2 = U. 

^From the above discussion, it is also easy to see that the interaction induced by the \Xij\ fluctuations only modify 



G e 



G 



Ui + 2i(cos k x + cos k y ) 
U 2 



(C14) 



In summary if we treat bosons as a free Bose gas and include attraction induced by gauge and \Xij\ fluctuations, 
the electron Green function is approximately given by Eq. C14. However, different treatments of fluctuations result 
in jj^ in a range from —1/3 to 1 and absolute magnitudes of U% t 2 are of order t. 
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-2 1 -2 1 

FIG. 7. The electron spectral function for U2/U1 = —1/3 and, from top down, (a) k = (— 7r/4, 7r/4) — » (tt/4, 3n/4), (b) 
fc = (-7r/8,7r/8) -> (3tt/8, 5tt/8), (c) fc = (0,0) -> (tt/2, tt/2), and (d) fc = (0,tt) -> (0,0). We have chosen J = 1. 




-2 1 -2 1 

FIG. 8. The electron spectral function for U2/U1 = 1 and, from top down (a) k = (— 7r/4, tt/4) — > (tt/4, 37t/4), (b) 
k = (-7T/8, tt/8) -v (371-/8, 5tt/8), (c) fc = (0, 0) -> (tt/2, tt/2), and (d) it = (0, tt) -> (0, 0). We have chosen J = 1. 

In Fig. Q and || we plot the electron spectral function for = —1/3, 1. We have chosen J = J/2, t = t = 2 J, 
X = 1, A/x = 0.4, x — 0.1, and T = 0.1J. The value U\ is determined from requiring the renormalized electron 
Green function to satisfies the sum rule 

"°du;d 2 k 

— —-ImG e = x (C15) 
The main purpose of the above study is to understand the ambiguity in the electron spectral function due to our 
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uncertainty in treating gauge fluctuation. We find that the gap near (0 ± n) and (±7r, 0) survives the inclusion of 
gauge fluctuation. However spectral functions near (±|,±|) are modified. For w < the electron spectral functions 
are quite similar for the three different choices of ^ = — 1/3, 0, 1. But for lj > the spectral functions show some 
notable differences. 
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